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Abstract 

^ . In this paper, we study the properties of baryons by using a holographic dual of 

QCD on the basis of the D4/D8-brane configuration, where baryons are described by 
a soliton. We first determine the asymptotic behavior of the soliton solution, which 
allows us to evaluate well-defined currents associated with the U(Nf) L x U(Nf) R chi- 
ral symmetry. Using the currents, we compute the static quantities of baryons such 
as charge radii and magnetic moments, and perform a quantitative comparison with 



experiments. It is emphasized that not only the nucleon but also excited baryons, 
such as A, iV(1440), iV(1535), etc., can be analyzed systematically in this model. We 
also investigate the form factors and find that our form factors agree well with the 
results that have been well established empirically. Using the form factors, the effec- 
tive baryon-baryon-meson cubic coupling constants among their infinite towers in the 
model can be determined. Some physical implications following from these results are 
discussed. 
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1. Introduction 



The gauge/string duality has opened up new technology to analyze strongly coupled 
gauge theories. Although it has not been proven yet, there is a lot of highly nontrivial 
evidence suggesting the duality between four-dimensional gauge theories and string theory 
in higher- dimensional curved backgrounds, at least for the cases with supersymmetry and 
conformal symmetry^"® (f or a review, see Ref. H])). It is obviously important to see if 
this idea can be applied to realistic QCD in many respects. One of the advantages of 
QCD is that we can use experimental data to check if the duality really works. In general, 
the analysis of non-supersymmetric strongly coupled gauge theories is very difficult, and 
hence it is almost impossible to establish such duality. However, for QCD, we can omit the 
complicated calculations on the gauge theory side and simply compare the results on the 
string theory side with the experimental data to obtain nontrivial evidence of the duality. 

If the gauge/string duality is really applicable to QCD, it will provide us with new deep 
insight into the theory of strong interaction in both conceptual and practical terms. It 
will tell us that both QCD and string theory in a higher-dimensional curved background 
(holographic QCD) can be a fundamental theory of strong interaction at the same time. Of 
course, in the high-energy weakly coupled regime, QCD is a better description. However, at 
least in the large-iV c strongly coupled regime, the string theory can be used as a powerful 
tool to calculate various physical quantities. 

In Refs. E]) andE]), it was proposed that QCD with Nj massless quarks is dual to type 
IIA string theory with iVj-probe D8-branes in the D4-brane background used in Ref. [7j) at 
low energy. In this model, mesons appear as open strings on the D8-branes and baryons 
are expressed as D4-branes wrapped on the nontrivial four-cycle of the background, as was 
anticipated in Refs. [S]) and[H]) in the context of the AdS/CFT correspondence. The effective 
action for the open strings turned out to be a five- dimensional U(Nf) Yang-Mills-Chern- 
Simons (YM-CS) theory in a curved background. Via the equivalence between the wrapped 
D4-branes and the instantons in the D8-brane world-volume gauge theory,^ the baryons 
were obtained as instanton configurations in the Euclidean four- dimensional space in this 
five- dimensional gauge theory^ (see also Ref. ITT]) ). 

The idea of realizing baryons as solitons was pioneered by Skyrme in the early 1960sP^ 
He started with an effective action for a pion with a fourth-derivative term (Skyrme term) 
and constructed a soliton solution (Skyrmion) in which the [7(iV/)-valued pion field carries 
a nontrivial winding number interpreted as the baryon number. In the model of holographic 
QCDP® this construction is naturally realized. It can be shown that the five- dimensional 
gauge theory is equivalent to a theory of mesons including pions as well as infinitely many 
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vector and axial-vector mesons. The pion sector is given by the Skyrme model and the 
instanton configuration reduces to a Skyrmion, realizing the old idea of Atiyah and Mantori^ 
who proposed expressing a Skyrmion using a gauge configuration with a non-zero instanton 
number. 

The Skyrme model was further developed by Adkins et al.,® who analyzed various static 
properties of the nucleon and A by quantizing the Skyrmion. It is natural to expect that the 
analysis can be improved by applying their strategy to the holographic model of QCD, since 
the model automatically contains the contribution from the vector and axial-vector mesons 
such as p, to, ai, etco In fact, the quantization of the instanton representing a baryon 
was studied in Refs. fTBjl - fTH]) . It was found in Ref. [TB]) that the spectrum of the baryons 
obtained in the model, including nucleons, A, iV(1440), iV(1535), etc, is qualitatively similar 
to that observed in nature, although the prediction of the mass differences tended to be 
larger than the observed differences. In Refs. [T7I) and fTHI) . a five-dimensional fermionic field 
was introduced into the five- dimensional gauge theory to incorporate the baryon degrees of 
freedom, where the interaction terms between the two fields were analyzed by quantizing 
the instanton. Using this five-dimensional system, the interaction of the baryons with the 
electromagnetic field and mesons was computed, which was found to be in good agreement 
with the experimental data. 

In this paper, we analyze the static properties of baryons, such as the magnetic mo- 
ments, charge radii, axial couplings, etc. of nucleons, A, iV(1440), and iV(1535), generalizing 
Refs. IT4"]) and [Ttjj) . These quantities are obtained by calculating the currents corresponding 
to the chiral symmetry that contain the information of the couplings to the external elec- 
tromagnetic gauge field. As we will see, most of our results are closer to the experimental 
values than the results given in Ref. IT4~j) for the Skyrme model. For some quantities of the 
excited baryons, our results are predictions. We also consider the form factors of the spin 
1/2 baryons and study how they behave as functions of the momentum transfer. 

Note that there are some closely related works that address the same subject. As men- 
tioned above, Refs. [T7|) and ITS) analyze the properties of baryons using the five-dimensional 
gauge theory with a five- dimensional fermion. On the other hand, our analysis is based on 
the quantum mechanics obtained through the usual procedure of quantizing solitons. These 
two approaches should in principle give the same results, and in fact we will find that some 
of our results are consistent with those in the studies, although the relation is not completely 
transparent. The work of Hata et al™ is more directly related to ours. The main difference 
is the definition of the currents. As discussed in their paper, their currents are problematic, 



*' See Ref. [T5]) for an attempt to include the effect of the p meson in the Skyrmion based on the 
holographic QCD. 
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since they are not gauge-invariant and have some ambiguities in the definition. 

The paper is organized as follows. In §2, after briefly reviewing the main points of the 
model and the construction of the baryon, we calculate the currents corresponding to the 
chiral symmetry. The applications are given in §3. We calculate various physical quantities 
including charge radii, magnetic moments, axial coupling, etc., and compare them with 
the experimental values for nucleons. We also present these values for excited baryons as 
our predictions. In §4, we compute the form factors of the spin 1/2 baryons and study 
their properties including cubic baryon-baryon-meson couplings. Section 5 is devoted to a 
summary and discussion with tables summarizing our results including our predictions for 
future experiments. Two appendices summarize some technical details. 



§2. Currents 



2.1. The model 



In this subsection, wegive a brief review of Refs. [5]), |6]), and fT6|l with an emphasis on the 
construction of baryonsJlj Based on the idea of gauge/string duality, it was proposed that 
the meson effective theory is given by a five- dimensional U (N/) YM-CS theory in a curved 
background. The action of the model is 



S — Sym + Scs 
Sym = — K I d^xdz tr ' 



h{z)J% + k{z)^ z 



24-k z ./ M 4 xI 



Here, v — 0,1,2,3 are four-dimensional Lorentz indices, and z is the coordinate of the 
fifth dimension. The quantity A = A a dx a = A^dx^ + A z dz (a = 0, 1, 2, 3, z) is the five- 
dimensional U(Nf) gauge field and J 7 = ^T^dx " A dx^ = dA + iA A A is its field strength. 
The constant k is related to the 't Hooft coupling A and the number of colors iV c as 

XN 

r = a\N c . (2-2) 



2167T 3 



Although it is not explicitly written in (12- 1[) . the mass scale of the model is given by a 
parameter Mkk, which is the only dimensionful parameter of the model. In Refs. |5]) and [6]), 
these two parameters are chosen as 

M KK = 949 MeV , k = 0.00745 , (2-3) 



*' Our notation is mostly consistent with Ref. [T5| except that we do not use the rescaled variables 
defined in (3.9) of Ref. US]). 
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to fit the experimental values of the p meson mass and the pion decay constant f w ~ 
92.4 MeV. In this paper, we mainly work with Mkk = 1 unit. The functions h(z) and 
k(z) are given by 

h(z) = (1 + z 2 )- 1 / 3 , k(z) = l + z\ (2-4) 
and ui 5 (A) is the CS 5-form defined asCl 

u 5 (A) = tr (aF 2 - l - A 3 J^-^A 5 ) . (2-5) 

The action (12 -ip is obtained in Ref. |5j as the effective action of iV/-probe D8-branes placed 
in the D4-brane background studied in Ref. [7]) and is thought to be an effective theory of 
mesons in four-dimensional (large N c ) QCD with Nf massless quarks. 

In this paper, we consider the Nf = 2 case, and the U{2) gauge field A is decomposed as 

^ = A + ^Y = Aa Y + A Y = ^ AC Y ' (2 ' 6) 

c=o 

where r a (a = 1, 2, 3) are Pauli matrices and r° = 1 2 is a unit matrix of size 2. Then, the 
equations of motion are 

Nr. 



- K (h{z)d v F^ + d z (k(z)F^)) + ^^ 2 "' a5 {Ka 3 K 4a5 + Fa 2 a 3 F a4as ) = , (27) 

- k (h(z)D v F^ + D z (k(z)Fn) a + ^ m2 '" at 'Ka 3 F<u a5 = , (2-8) 

- k k{z)d v F™ + j^e^ 5 [f;^F^ + F, 2 , S F^ 5 ) = , (2-9) 

Nc 

64vr 2 ' 



-Kk{z) (D u F zv ) a + ^""V^ = , (2-10) 



where D a = d a + iA a is the covariant derivative. The baryon in this model corresponds to a 
soliton with a nontrivial instanton number on the four- dimensional space parameterized by 
x M (M = 1, 2, 3, z). The instanton number is interpreted as the baryon number Nb, where 

NB = 64^ j ^ xdz e M 1 M 2 M 3 M 4 FM 1 M 2 F M 3 M 4 ■ ( 2-11 ) 

Unfortunately, because the equations of motion are complicated nonlinear differential 
equations in a curved space-time, it is difficult to find an analytic solution corresponding to 
the baryons. However, as observed in Refs. IToT) and fTTl) . the center of the instanton solution 
is located at z = and its size is of order A~ 1//2 , and hence we can focus on a tiny region 



*' Here we omit the symbol 'A' for the wedge products of A and T (e.g. AT 2 — A A J- A J-) in the CS 
5-form. 
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around z = for large A, in which the warp factors h(z) and k(z) can be approximated by 1. 
Then it follows that the static baryon configuration is given by the BPST instanton solution 
with the £7(1) electric field of the form 



Ad 



^0 



N c 1 

87T 2 K£ 2 



1 - 



(p 2 + e 



2^2 



An = A 



M 







(2-12) 



Here 



/(0 



e+P 2 



g{x) 



z) 



i\x 



X)-T 



zy 



\x-X\ 



(2-13) 



with X M = (X 1 , Jf 2 , Jf 3 ,Z) = (X,Z) being the position of the soliton in the spatial M 4 
direction and the instanton size p. Substituting this configuration into the action and taking 
the nontrivial z dependence of the background into account as a 1/A correction, we find that 
p and Z have a potential of the form 

- 2 iV 2 1 



U(p, Z) 



lli^K 1 



6 5(8vr 2 fi;)V 



Z 2 
3 



(2-14) 



which is minimized at 



Pel 



8ti 2 k 



. 



(2-15) 



In order to quantize the soliton, we use the moduli space approximation 
methocP^'''^ and then the system is reduced to the quantum mechanics on the instanton 
moduli space. The 577(2) one-instanton moduli space is simply given by Ai ~ M 4 x M 4 /Z 2 , 
which is parameterized by (X, Z) and yi (I = 1, 2, 3, 4) with the Z 2 action yj — > —yi. The 
size of the instanton p is related to yi by p = \/ y\ + • • • + yi, and aj = yi/p represent the 
SU(2) orientations of the instanton. The Lagrangian of the collective motion of the soliton 
was obtained in Ref. [16]) as 



M 



(X 2 + Z 2 ) + M yj-U(p, Z) 



(2-16) 



where M = 8tt 2 k. The Hamiltonian is given by 



H 



with the canonical momenta 



1 



2M 



(P 2 + P|) 



AM, 



■n] + u( P ,z) 



M X 



._d_ 
'dX 7 



MnZ 



._d_ 

'dZ 



/J 7 = 2M yi 



._d_ 



(2-17) 



(2-18) 
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This system is manifestly invariant under 5*0(4) rotation acting on y 1 . As argued in 
Refs. E]) and H~6"|) the generators of SO (A) ~ (SU(2)j x S77(2)j)/Z 2 symmetry correspond 
to the isospin and spin operators given by 

% ( d d d \ i ( d d d \ 

2 V <% a 02/4 <9y c / 2 \ dy a dy 4 dyj 

(2-19) 

respectively 

The explicit eigenfunctions of the Hamiltonian (12 • 1 7|) are obtained in Ref. fl~6|) . They are 
characterized by the quantum numbers (I, I3, s,n p ,n z ) as well as the momentum p. Here 
/ = 1, 3, 5, ■ ■ • are positive odd integers related to isospin / and spin J by I = J = 1/2. Note 
that (12- 19|) implies I 2 = J 2 , and hence only the states with I = J appear in the spectrum as 
in Ref. H4"]) for the Skyrme model. ^3 and s denote the eigenvalues of the third component of 
the isospin and spin, respectively. n p and n z are nonnegative integers corresponding to the 
excitations with respect to p and Z, respectively. The species of the baryon are specified by 
B = (I, I3, n p , n z ) and the wavefunctions are of the form 

\p,B,s) = \p)\B,s) , (2-20) 
with \p) = ( 27 r)3/2 etp ' X " F° r l a ^ er use 5 we normalize the baryon state as 

(B,s\B',s'} = 5 B>B ,5 s>s , . (2-21) 

For example, the proton and neutron are interpreted as the particles with B = (1, 1/2, 0, 0) 
and B = (1, —1/2, 0, 0), respectively, and the corresponding wavefunctions with s — 1/2 are 

\pt)(xR(p)^ z (Z)(a l + ia 2 ) , \n f) oc R{p)^ z {Z){a A + ia 3 ) , (2-22) 

respectively, where 

R(p) = p -^V^^e-^ p2 , MZ) = e~^ z2 , (2-23) 

up to normalization constants. See Appendix IA. II for more details. 

The excitation numbers (n p , n z ) are the quantum numbers which Skyrmions cannot have, 
and are thus peculiar to the instanton picture of baryons obtained in Ref. [16]). For example, 
B = (1, ±1/2, 1,0) corresponds to the Roper excitation JV(1440), and B = (1, ±1/2, 0,1) 
corresponds to iV(1535). Higher spin baryons are also included, such as B = (3, 1 3 , 0, 0) with 
h = ±3/2, ±1/2 giving Z\(1232). 
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2.2. Currents 

As in the usual effective field theory approach to the hadrons, it is useful to introduce the 
external gauge fields Al, and Ar, associated with the chiral symmetry U(Nj)l x U(Nf)n- 
The currents and associated with the chiral symmetry are then read off from the 
terms linear with respect to the external gauge fields in the effective action as 

A R ,J») . (2-24) 
J£ , (2-25) 



V<->J2) , (2-26) 

(2-27) 

are external vector and axial-vector fields. 

We can apply this idea to compute the currents in the present model, which leads us to 
a prescription analogous to that used in the AdS/CFT correspondence as first explored 
in Ref. [11]) and discussed recently in Ref. [19]). As studied in Refs. [5) and [6]), the chiral 
symmetry is identified with the constant gauge transformation at infinity z = ±oo, and the 
external gauge fields Al, and Ar, can be introduced by considering the five-dimensional 
gauge field with the boundary conditions 

A^af, z -> +oo) = Al^x") , A^, z -> -oo) = A Rfl (x») , (2-28) 

respectively. To compute the currents, we substitute a solution of the equations of motion 
with this boundary condition into the action and keep only the linear terms in the external 
gauge fields Al and Ar assuming that the external fields are infinitesimal. Consider the 
gauge configuration of the form 

A*(af, z) = At(x», z) + 6Aa(x^ z) , (2-29) 

where A& is a classical solution of the equations of motion with A^(z = ±oo) = and 
SA a {x, z) is an infinitesimal deviation from it satisfying 

SA^x", z -> +oo) = A LfM (x") , SA^af, z ->■ -oo) = Ar^) . (2-30) 



S\ 0{ A L ,A R) = - 2 J d ^ { A ^l 



In terms of the vector and axial- vector current 

Uy — U L T Or , U A — U L 

f[2 T 24l) becomes 



S\ 0{AL>AR) = -2jd*xtr (vWj; 

where 

} = \ (A L , ± A R ,) 



8 



Substituting this configuration into the action (12-11) . we obtain 



S \o(A L ,A R ) - K I rf4x2tr 



tr (5A»k{z)rf z ) 



J z=— oo 



which implies 



J L , = -K (k( Z ) J&) 



-K 



k(z)T. 



z=+oo 
-\ z=+oo 



Jrv = +« (k(z) rf z ) 



An 



-K 



%{z)k(z)F 



z=— oo 

2 = + 00 

cl 

fXZ 



(2-31) 

(2-32) 
(2-33) 



where ipo( z ) — ~ arctanz. Here only the surface terms at z — ±oo remain in (l2-3ip because 
of the equations of motion. 

Note that the currents in (I2-32P are consistent with the four- dimensional effective action 
obtained in Ref. |6). It was shown in Ref. |6])|*^ that the four-dimensional effective action 
derived from our model has the following terms: 



S 



0(Al,Ar) 



where H (x), v™(x), and a™ (a?) are the pion, vector meson, and axial- vector meson fields, 
respectively, and /„., g v n, and g a n are the decay constants of these mesons, respectively. 
These mesons are related to the five-dimensional gauge field as 



A„(x,z) = ^t»;(x)^2n-i(^) + J2 a ^ x )^(z) , A z (x,z) = n{x)<j) {z) , (2-35) 



n=l 



n=l 



where 0o 



1 1 



//ot k(z) 

the eigenfunctions of the eigenequation 



and {ip n {z)} n =i,2,- is a complete set of the functions of z consisting of 



h(z) 1 d z (k(z)d z i[) n ) = A„Vv, 



with the normalization condition 



k I dzh{z)ip n ip m = 5 r 



(2-36) 



(2-37) 



We can show that ip2n-i(z) an d ^2n{z) are even and odd functions of z, respectively, if we 
arrange ip n (z) such that the eigenvalues satisfy Ai < A 2 < A 3 < ■ • ■ . The decay constants 
are given by 



Z = + 00 



g an = -2k {k(z)d z i) 2 n) 



I z=+oo 



(2-38) 



*) See Eq. (5.12) in Ref.©. 
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Note that the functions if> n (z) behave as 0(z 1 ) in the z — > ±00 limit and the decay constants 
g v n and g a n are determined as the coefficients in front of 1/z as 

1p2n-l{z -> ±00) ~ , ^2n(^ ~> ±Oo) ~ |^ . (2-39) 

ZKZ zkz 

The following expressions are also useful: 

9v" = A 2n -i« / dz h(z)ip 2 n-i , 9a« = A 2ri K / dz h(z)tfj 2n ipo , (2-40) 



which are obtained by using (I2-38P and (l2-36p . 

Comparing (I2-26P and (l2-34j) . the vector and axial- vector currents are obtained as 



n=l n=l 

Using (I2-35P and (12 -38 p . we can easily check that f!2-33|) and (!2-4ip are equivalent. Note that 
the vector current is expressed by the vector mesons v ™, which is a direct consequence of the 
complete vector meson dominance of the model found in Ref. |6|). (See also Refs. [TT]) and 



As in f)2-6p . we decompose the currents as 

rC 

T 



c=o 

Then the baryon number current is given by 



2 ~ 2 



k(z)F" z . (2-43) 

As a check, the baryon number density is computed as 



= "oi / dze ° MlM2M3M4F MiM 2 F a M 3 M 4 + (total derivative) , (2-44) 
where we have used the equation of motion (12- 7p . This is consistent with (12-lip as expected. 

2.3. Asymptotic solution 

In order to calculate the currents (12-321) . we have to know how the field strength T^ z 
behaves at z = ±oo. Unfortunately, we cannot directly use the solution (I2-I2p since it is 
only valid in the region £ <C lo In this subsection, we study how to extend the solution to 
the 1 <^ £ region. 



*' Here we are assuming Z ~ Z c \ = 0, since the expectation value of f(Z) is approximated by its 
classical value f{Z c \) for the large N c and large A limit. 
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Let us first summarize the gauge configuration in the ^ < 1 region. Here we include the 
time-dependent moduli parameters, which are treated as operators in quantum mechanics 
(12- 17j) . As explained in Ref. [16]), the SU(2) gauge field takes the form 

A M = VAfjV- 1 - iVduV- 1 , (2-45) 

where A M is the solution ( I2-I2p and V satisfies 

- iv- l v = -x m a m + x a f(0 g^ 1 , (2-46) 

which follows from the Gauss law constraint. Here 

X a = -i tr (r a a-'a) = -± tr (r^y) = -J— 2 J a , (2-47) 

where y = y 4 + iy a r a , a = a 4 + ia a r a = y/p, and J a is the spin operator defined in (I2-19I) . 
It is convenient to perform the gauge transformation 

A a —> A^ = GA a G- 1 - iG8 a G~ l , (2-48) 

with G = ag~ l V~ l . Then, 

= - f(0) a a 1 + i(l - /(£)) X M a (g-'dMg) a 1 , (2-49) 
A% = -z(l - /(0) a (g-'dMg) a 1 . (2-50) 

This choice of gauge is useful in considering the asymptotic behavior, as we will see in the 
following, while a singularity develops at £ = 0. 

The U(l) part is treated as a perturbation in the background given by the SU{2) gauge 
configuration obtained above. The leading contribution to the U(i) part is obtained by 
solving the following linearized equations of motion (in the Lorenz gauge): 

9mO m A° = 9 3 : ^ ^ , (2-51) 



ir 2 a\ (£ 2 + p 



jM ?i _ 3 p f vi , X /Jai„i xia„\ , P X 



B^X = ^ + ^,r> - rz) + (2-52) 

OmVZ = 4t77^7 Z + **f + ^). (2-53) 



7r 2 aA (£ 2 + p 2 ) 4 V 2 p 

Here we substituted ( I2-49P and ( l2-5(jp into the equations of motion (12 -7p and f!2-9p with the 
warp factors h(z) and k(z) approximated by 1. We neglect the terms including 8q, because 
we are interested in slowly moving solitonso 

*' The time derivative squared such as Z and Z 2 can be traded by Z or Z 2 via its Schrodinger equation 
or the relation Z — Pz/Mq; thus for large N c and A, it is approximated by its classical value, which vanishes. 
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The regular solution is found to be 



A 
A 

A, 



8n 2 a\ e 
1 



P 1 



iir 2 a\ 
1 



e + 2 P 2 



i e + 2 P 2 

8vr 2 aA (£ 2 + p 2 ) 2 

2 



]ir 2 a\ 



(e + p 2 ) 2 
e + 2 P 2 



x l + 



z + 



p 



{e + p2)2 

,2 



X 



(e m3 x J 



5 m z) + 



px 
P 



P 



X a x a pz 
~^ + ~P 



(2-54) 
(2-55) 
(2-56) 



ue+p 2 ) 2 (e+p 2 ) 2 

Note that (I2-55P and (I2-56P are neglected in Ref. ITB"]) . since the energy contributions from 
them are subleading in the 1/A expansion given in Ref. 16). However, here we keep them 
because they give the leading contribution to the isoscalar current density. It can also be 
shown that the F 2 terms in the equations of motion (I2-7|) and (12 • 9 j) are subleading in the 
1/A expansion, justifying the above perturbative treatment for the U(l) part. 

So far, we have established a solution that is valid in the region £ <C 1. We now consider 
how to find the solution in the 1 <C £ region. The key observation is that all the components 
of the gauge field in (1249]) . (E2D}, (EMJ), (E3S}, and fEFBT)]) are suppressed in the p < £ < 1 
region in the large A limit. This implies that the nonlinear terms in the equations of motion 
can be neglected in this region for large A. Our strategy is to find a solution of the linearized 
equations of motion in the p <C £ region that smoothly connects the previous solution in the 
overlapping region p <C £ < 1. 

For this purpose, we note that for p < £ < 1, the gauge field ( !2-49p . fl2-50p . f l2-54p . 
f!2-55p . and ( I2-56P are approximated as 
1 



A, 
M 

A, 



2a\ 



G aat (x,z;X,Z) 



2a\ 
1 

2oA 



Z + 



r 

2 

p 2 ( x a d 



_iaj 







+ 



dXi 
p d 



dZ 



+ 



p d 



pdX i 



G aat (x,z;X,Z) 



2 V 2 dX a p dZ 



G aat (x,z;X,Z) 



(2-57) 
(2-58) 

(2-59) 



Aq ~ 4vr 2 p 2 ia a" 1 G flat (f, z; X, Z) 



2ir 2 p 2 aT a a- 1 



X 1 



d 



-iaj 



dXi 



■ d 

rat 

6 dZ 



d 



dX a 



G flat (f, z;X,Z) 



A" 



G 



2tt 2 p 2 



ara — + e iia ar a — 
~ -2n 2 p 2 aT a a- 1 ^—G R£lt (x,z;X,Z) , 



G flat (f,z;X,Z) 



where 



G nat (x,z;X,Z) 



1 1 



(2-60) 
(2-61) 
(2-62) 

(2-63) 
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is the Green's function in the flat IR 4 that satisfies 

d M d M G^\x, z; X, Z) = 5 3 {x - X)5{z - Z) . (2-64) 

We can easily check that the gauge configuration (1 2 • 5 7 1) (1 2 • 6 2 p satisfies the Maxwell and the 
linearized YM equations without sources: 

as well as the gauge condition 

d a A a = , d a A G = . (2-66) 

In order to connect this solution to the large £ region, we have to take into account the 
effect of the curved background. The linearized equations of motion and the gauge condition 
that generalize ( I2-65P and (12-66P to the case with nontrivial warp factors h(z) and k(z) are 

h(z)%Ai + d z (k(z)d z Ai) = , d% + d z (h{z)- x d z {k{z)A z )) = , (2-67) 

h(z)dlAf + d z (k(z)d z Af) = , d\A G z + d z (h(z)- l d z (k(z)A G z )) = , (2-68) 

and 

h{z)d^A fX + d z (k(z)A z ) = , h(z)d»A G + d z (k(z)A G ) = , (2-69) 

respectively. 

To solve these equations, we define Green's functions in the curved space as 

oo 

G(x,z;X,Z) = kJ2Mz)MZ)Yu(\x-X\) , (2-70) 
n=i 

oo 

H(x,z;X,Z) = K J2<f>n(z)<l>n{Z)Yn(\Z - X\) , (2-71) 

n=0 

where {i>n(z)}n=i,2,- is the complete set defined in (I2-36P and (12-371) . {4> n (z)}n=o,i,- is an_ 
other complete set given by 

o (z) = _L_L <f> n {z) = ^=d z il> n {z) , (n = l,2,...) (2-72) 
and Y n ( r) is the Yukawa potential with meson mass fn n — \J A n , 

y - (r) = "s— • (2 ' 73) 
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which satisfies 

- X n )Y n (\x - X\) = 5 3 (x - X) . 
Note that the normalization of the functions {</> n }n=o,i,2,- is fixed by 

« y dzk{z)(j) n (j) m = b mn ■ 

Using (12 -36 p . (!2-72j) . (12- 74j) . and the completeness conditions 

oo oo 

Kh(z)^2^ n (z)ijj n (Z) = 5{z - Z) , Kk{z)y^ j (j) n {z)(j) n (Z) = 5{z - Z) 

n=l n=l 

it is easy to verify 

h{z)d 2 G + d z {k{z)d z G) =5 3 {x- X)5{z - Z) , 
d 2 H + d z {h{z)- 1 d z {k{z)H)) = kiz)- 1 5 3 {x - X)8(z - Z) , 
d z (k(z)H) + h(z)d z G =0 . 



(2-74) 



(2-75) 



(2-76) 



(2-77) 
(2-78) 
(2-79) 



Then, the solution of equations ( l2-67[) ( l2-69[) is obtained by replacing the Green's function 



G flat in (ESZD-dZiSI with G or H as follows: 



A 



2a\ 



G(x, z; X, Z) , 



1 



2a\ 
1 

2^A 



X l + 



p- X 



MJ 







dXi 



t 9_ 
dZ 



+ 



p d 



pdX i 



G(x,z;X,Z) 



v + ff_( X^_d_ pd_ 
2 V 2 dX a pdZ 



H(x,z;X,Z) 



(2-80) 
(2-81) 

(2-82) 



Aq ~ A-K 2 p 2 ia a 1 G{x, z\ X, Z) 



dX a 



G(x,z;X,Z) 



A 



G 



0^2 Jl-i^a _ — 1 

-Z7T p ar a 



A 



G 



d 

-27T 2 p 2 aT a a- 1 -—H(x,z]X,Z) . 

f Q X a V » J 



(2-83) 
(2-84) 
(2-85) 



Here, we neglected the terms including &q as before. 

Since the Green's functions G and H approach G flat for /) < ( < 1, this solution is 
smoothly connected with the previous solution (12-57^ (12-62^ in this region, as expected. 
Note also that the Green's functions G and H vanish at £ — > oo, and hence the linear 
approximation of the equations of motion does not break down all the way to infinity. 
Therefore, we can read off the behavior of the gauge field at z — > ±oo, which is needed to 
calculate the currents, from (|2-80p - (j2-85p . 
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2.4. Computation of the current 

Since the wavefunctions of the low-lying baryon states, such as (j2-22j) . are dominant for 
Z ~ 0(X~ 1 ^ 2 N c 1 ^ 2 ) <C 1, we can use the approximation h(Z) ~ k(Z) ~ 1 and the relation 

- d\i) n {z) ~ \ n Mz) , 

which follows from (!2-36[) . This implies 

8 Z H + d z G ~ , + <9f )G ~ , + df)# ~ . 

Using the asymptotic solution ( I2-8Q[) (!2-85[) and the relations f l2-87p . we obtain 
1 



F 0z 
F 



2a\ 
1 

2aX 



Z^tf - X%G - ^j- ((9i9 a - 5 ia <9|)# - e^djd.G) 



F 0z ~ 2n 2 d {p 2 aT a a- 1 )d a H - An 2 p 2 iaa' x d z G 

- 27r 2 p 2 aT a a~ 1 X i {{d a di - 5 ia d 2 )H - e^d 3 d z G) 
F iz ~ 27r 2 p 2 ar"a- 1 ((5^ a - <f a ,9j)# - e^d&G) 

for Z < 1 < 2. 

The following formulas are also useful: 



G v (Z,r) = k(z)d z G 
G A (Z,r)= [^{z)k{z)c>G 
H v (Z,r) = 



J z=— oo 



= - ^ 9v^2n-\ {Z)Y 2n _ x (r) , 

n=l 

oo 

;=+oo - 

= - > g a ^2n(Z)Y 2n (r) 

:=— oo 

n=l 

oo 



n=l 



n-H' ; , 



H A (Z,r) 



%{z)k{z)H 



z=+oo 



1 1 1 



J z=— oo 



27T 2 fc(Z) r ^ A 2 n 
n=l 



2n r 



where r = |af — X|. Note that and H A are even functions with respect to Z, while G A 
and if y are odd functions. 

Now, we are ready to write down the currents from (12-321) and (!2-33[) . The vector and 
axial-vector currents are obtained as 



tD _ N c V A 
j v,a — ~~2~ U ' 



V,A 



iVc 

2 



ZdiH v ' A - X l G VA - P —f- ((did a - 5 ia d])H v ' A - e iaj djG v ' A ) 



(2-96) 
(2-97) 
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J° A = 2tc 2 k d {p 2 aT a a- l )d a H v > A - 2p 2 ia a,- l G v > A 

- fa^a^X* ({dad, - 6 ia cP)H v > A - e iaj djG v > A ) 1 , (2-98) 
r vA = -2n 2 K p 2 ^^- 1 ({did a - 5 m d 2 )H V ' A - e iaj d 3 G v ' A ) . (2-99) 

Using the relation ()2-87p . it is easy to check that these currents are conserved (up to terms 
including <9q)l2 

§3. Static properties of baryons 

In this section, we study the static properties of baryons as applications of the currents 
obtained in the previous section. 

3.1. Baryon number density, isoscalar mean square radius 

The baryon number density is obtained from (I2-43|) and (I2-96P as 

oo 

J B = -X>"^2n-l(Z)y 2n _i(r) . (3-1) 
n=l 

As a check, the baryon number charge is calculated using this expression as 

poo 

N B = dr Anr 2 (j° B (r)) (3-2) 



o 

oo 

n=l A2 ™- 1 

= 1 , (3-3) 

where (O) = (B, s\0\B, s) is the expectation value with respect to a baryon state \B,s). 
Here we have used (l2lH and fl2 7 76|) . 

The baryon number density per unit r is given by the integrand of (13- 2p . 



p B (r) = Anr 2 (j°(r)) = r (^n-i(Z)) e~V^ T . (34) 

n=l 

Then the isoscalar mean square radius is 

oo 

« — , n „ 

(lp2n-l{Z)) • (3-5) 



(r 2 ) I=0 = drr 2 p B (r) = 6J2^T— 

JO n=1 A 2n-1 



In the large N c and large A limit, the baryon wavefunction is localized at Z = 0, as 
we can see from the wavefunction f)2-22p for the nucleon states, and hence the expectation 



Here we neglect the effect of the U(1)a anomaly, since it is subleading in the 1/N C expansion. 
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value (ip2n-i(Z)) can be approximated by its classical value ^% n -\{Z c \) = ^ 2 n-i(0). In this 
approximation, it is possible to evaluate the isoscalar mean square radius (I3-5P as follows. 
Note that the function 



n=l 



A 



2n-l 



satisfies 

- d z (k(z)d z F(z)) = 6 h(z) , F(z) = F(-z) (3-7) 

and the boundary condition F(z) — > {z — > ±oo). The first relation is obtained from (!2-36j) . 
(T23QD, and fP2~76]) . The solution of (EI} is given by 

F{z) = F - [ dz'kiz')- 1 [ dz"6h(z"). (3-8) 
Jo Jo 

The constant Fq is fixed by the boundary condition and we obtain 

F = dz'k(z')- 1 dz"dh(z") ~ 14.3 . (3-9) 
Jo Jo 

Therefore f!3-5p can be evaluated as 

( r2 ) I= o = ( F ( Z )) - = F ° - 14 - 3 / M kk , (3-10) 
and if we use the value ( 12 -3p for Mkk, we obtain 

( r2 )J=o - °- 785 fm • (3' 11 ) 

The experimental value is (r 2 ) 1 /^ | ex p — 0.806 fmp3 and the prediction of the Skyrme model 
in Ref. Q3J is 0.59 fm. 

It is interesting to note that this value is independent of A and N c . The iV c independence 
is consistent with the analysis of baryons in large N c QCD.^ 1 The A independence suggests 
that the size of the baryon number distribution is governed by the scale of the vector meson 
mass rather than the size of the soliton p c \ ~ (9(A -1 / 2 ) for large A. 

Given the wavefunction of the baryon state, it is also possible to evaluate the expectation 
value (F(Z)) numerically!*)] For the nucleon wavefunction given by fl2-22j) . we obtain 

( r2 )J=o - °- 742 fm • (3" 12 ) 

*' Note that in the previous section the chiral currents were obtained using the approximation Z ~ 0, 
and thus the classical value for Z was used. Here, we simply assume the same chiral current and evaluate it 
using the quantum states of the baryons. 
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This value is the same for the other states with n z = 0, such as Z\(1232) and iV(1440), since 
Z dependence of the wavefunction is the same for these states. For the states with n z = 1, 
such as iV(1535), we obtain 

(r 2 ))^ ~ 0.699 fm . (3-13) 

See Appendix IA.2I for more details. 

3.2. Isoscalar magnetic moment 

The isoscalar magnetic moment is defined as 

/4 =0 = ^ e ijk J d 3 xx j J% = ^e ijk J d 3 xx j J^ . (3-14) 

It is easy to see that only the last term in ( I2-97P contributes to the integral. Then by 
evaluating the angular integral in (13- 14j) we obtain 

/4=o = J™dr4nr%J° B (r) . (3-15) 
Using f)2-47p . this can be evaluated as 

= £>L = J— (346) 

4 2M v ' 

where J 1 is the spin operator (12-19)1 . For example, for nucleon states with up spin, this gives 

( Pt I /4=o I Pt> = ( nt I /4=o I nt) = ^S 3i . (3-17) 



The g factor is defined as 



where a 1 is the Pauli matrix that acts on a spin doublet, and is the nucleon mass. If we 
use ( 12-3)1 and the experimental value for the nucleon mass ~ 940 MeV, we havfl 

3/=o = 9 P + 9n = M N /M ~ 1.68 . (3-19) 

The experimental value is gi = o | oxp ~ 1.76^ and the prediction obtained from the Skyrme 
model in Ref. 114)1 is 5 | /=o|anw — 



*> Note that is not calculated in Ref. [T5| because the total contribution from the zero point energy 
of the fluctuations around the soliton solution is difficult to evaluate. Here the g factor is computed simply 
to express the magnetic moments in the unit of l/(4Mjv)| exp . 

**) The same expressions as (|3- 171) and (|3-19[1 were obtained using the current in Ref. . However, the 
numerical values were different, because the values of Mkk and k used in Ref. HU) were different from ours. 
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3.3. Isovector charge density and charge radii 

The isovector charge density is given by (I2-98P and the isovector charge is evaluated using 
O as 

Qv = J d 3 xJ^ = -Air 2 Kp 2 iaar l . (3-20) 

A relation similar to ( I2-47|) . 

- i tr (r"a a" 1 ) = -± tr (r°y yt) = _L_/ a , (3 . 2 1) 
where I a is the isospin operator defined in (12- 19j) . implies 

Qv = Ia, (3-22) 

as expected. 

The isovector charge density per unit r is proportional to the angular integral of the 
isovector charge density with the normalization condition J °° dr p/=i(r) = 1. This turns out 
to be identical to the baryon number density: 

P/=i(r)=/> fl (r) . (3-23) 

Therefore, the isovector mean square charge radius (r 2 ) J=1 = J °° r 2 pj = i(r)dr is the same as 
the isoscalar mean square radius (r 2 ) I=0 evaluated in §3.11 This result is somewhat puzzling 
since it is known that the isovector mean square radius is divergent in the chiral limit 
There is however no contradiction. This divergence is due to the IR divergence of pion 
loops.®' Our analysis only involves a string world-sheet with disk topology, and hence the 
pion loops are not included. Therefore it will be important to include the quark mass in the 
models and estimate the contribution from the annulus diagram to make a comparison with 
the experimental value, which is beyond the scope of the present paper. 
The electric charge is defined by 

Qem = h + ~f , (3-24) 

which gives Q em = 1 for a proton (I 3 = 1/2) and Q em = for a neutron (J 3 = —1/2). Then, 
because of identity (13- 231) . the electric charge density is given by pe = Pi=o for a proton and 
Pe = for a neutron, and we obtain 

(r 2 ^ Sp = (r 2 ) / _ (for a proton) , (r 2 ^ En = (for a neutron) . (3-25) 



*' See Refs. I?7 [) -IM 1) for recent developments toward the incorporation of the quark mass in the model. 
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The experimental values in Ref. I23j) are 

Although our calculation does not reproduce the experimental value of the electric charge 
radius of the neutron (r 2 ) En , the vanishing of the neutron electric charge density seems to 
be a good approximation for reproducing the observed behavior of the electric form factor, 
as we will study in §4. 

The charge radius for the excited baryons is similarly found to be the same as their 
isoscalar mean square radius (r 2 ) I=0 obtained in §3.11 In particular, our analysis predicts 
that the Roper excitation iV(1440) has a charge radius equal to that of a proton (13- 12|) . while 
that of JV(1535), (EFTS]) , is smaller. 

3.4. Magnetic moment 

The isovector magnetic moment is defined as 

// /=1 = ^ e ijk J d 3 xx j tr( j£r 3 ) x 2 . (3-27) 

Here the additional factor of 2 in the integrand is due to our normalization of the current. 
Substituting (I2-99P into this expression gives 

/4 =1 = -47T 2 Kp 2 tr(ar l a" 1 r 3 ) . (3-28) 

For the baryon states of I — J — 1/2, we can use the identity^ 1 

(B\s'\tv(aT l a- l T a )\B,s) = —(a*) a , s (r a ) ^ , (3-29) 

with a 1 and r a being the Pauli matrices corresponding to spin and isospin, respectively. Here 
we have used the notation 



(o*)y. = x\ sl) ° % X(s) , (r c )i> h = r c Via , (3-30) 
where X(s) an d i>i 3 are defined as 

X(i/2) = ^3=1/2 = ( J J , X(-i/2) = tpi 3 =-i/2 = ( J J • (3-31) 
Then, we obtain 

(Pt I /4=i I Pt> = -( nt I /4=i I nt> = ^ </> *° , (3-32) 
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where (p 2 ) is the expectation value of p 2 for the nucleoli states. Using the wavefunction 
(HH2D, (p 2 ) is obtained as 



jdpp 5 R(p) 2 v ^+2 v ^nv? 2 

< P > = jdpp*R(p) 2 = 2N C '<* ■ (3 - 33) 

In the large N c limit, it can be approximated by its classical value p\ in (I2-I5p . However, for 
N c = 3, (13-331) implies (p 2 ) np=0 — 1.62 p^, suggesting that the l/N c corrections are relatively 
large for this quantity. 

Consequently, we obtain 



<Ptl/4=ilPT> = -Ml/^lnt) = l + 2 ^ N2J ^ ■ (3-34) 

Here we have recovered the Mkk dependence by dimensional analysis. The isovector g factor, 
defined as in (13-181) . is then 



gi=i =g P -g n = tt ^ 7.03 (3-35) 

for N c = 3, M N ~ 940 MeV, and fl23J. The experimental value is #/=i| eX p ^ 9.41 23 ' and the 
prediction obtained from the Skyrme model in Ref.[H]) is <?;=i|anw — 6.38. If we approximate 
(p 2 ) by its classical value p 2 1; we have gi=\ — 4.34, which is considerably smaller than the 
experimental value. It is interesting to note that going to (p 2 ) from p\ by multiplying by 
the ratio (iV c + \/5/2)/N c for large N c , as can be seen from ()3-33p . has an effect similar to 
the U N C — > N c + 2" rule discussed in Ref . ITT]) . The value of the anomalous magnetic moment 
obtained using this rule from the five- dimensional effective spinor field theory approach in 
Ref. [T7]) is close to oursO 

The magnetic moments for the proton and neutron (measured in units of the Bohr mag- 
neton /ijv = l/(2Mjv)) are given as 

P P = y = J\9i=o + 9i=i) , fi>n = y = j{9i=o ~ 9i=i) • (3-36) 

If we insert the values C I3-19[) and (13 -35 p . we get 

p p ~ 2.18 , p n ~ -1.34 , (3-37) 

while the experimental values are p p | ex p — 2.79 and p n | ex p — —1.91. Note, however, that 
since the N c dependences of the g factors and the magnetic moments are 

g I=0 ~ O(l) , g I=1 ~ 0(iV 2 ) , p /=0 ~ 0(1/JV C ) , p /=1 ~ 0(N C ) , (3-38) 



*) Again, the same expressions (|3-32[) and (I3-35|) can be found in Ref. but with a different current. 
Eq. (I3-32|) also agrees with zi^ an in Ref. 17), if we use the classical value (|2-15p for (p 2 )- 



21 



the contribution of the isoscalar component will be buried in the 1/N C corrections of the 
isovector component in linear combinations such as (13 -36 p . Therefore, it is more meaningful 
to consider gj =0 and gi=\ rather than g p and g n in our analysis. 

Let us consider excited baryons. For spin 1/2 excitations of the baryons, the magnetic 
moment remains the same as that of the proton/neutron if n p = (for example, iV(1535)), 
because it gives the same (p 2 )- On the other hand, the Roper excitation iV(1440) has n p — 1 
for which, with N c = 3, 



(3-39) 



See Appendix IA.2I for the details. Substituting it into f)3-32p and combining the result with 
the value of pi=o (which is the same value as that of the proton/neutron), we obtain, for the 
Roper excitation, 



p v * ~ 2.99 , p n * ~ -2.15 



(3-40) 



measured in units of the Bohr magneton of the nucleon l/(2Mjy). Our model predicts that 
the magnetic moment of the Roper is larger than that of the proton/neutron. 

For spin 3/2 baryons such as A, we need to reevaluate the matrix elements of the spin 
operator J\ (I3-29P and (I3 33[) . We find (details are described in Appendix I A. 3ft 



ti = o(B, h = 3/2) = 3/4 =0 (p t) , (B = A ++ , A + , A , ZT) (3-41) 

p) =1 (A ++ , J 3 = 3/2) = j?U(p t) , ^/=i(^ + , ^3 = 3/2) = t) , 

3r 9r 
t4 =1 (A°, J 3 = 3/2) = --/4 =1 (p t) , A«j=i(^". ^ = 3/2) = - T /iLi(P t) , (3-42) 



where c = (p 2 ) l=3 / (p 2 ) l=1 ^ 1.34 



Therefore we obtain 



p A ++ ~ 5.50 , pa+ ~ 2.67 , pao ~ —0.15 , p^- 



-2.97 



(3-43) 



using (I3-19P and f l3-35p . The experimental values for A ++ and A + ar 



PA++ |exp - 3.7 - 7.5 , p A A 



2.7t\i± 1-5 ±3 



|cxp 



(3-44) 



with which our result fl3-43[) is found to be consistent. A recent lattice result gives^ 



| lattice - 4.99 , p A + | lattice ~ 2.49 , p A o 



lattice - 0.06 , p A - (lattice — 



-2.45 , (3-45) 



*^ This factor was missing in the earlier versions of the present paper. We thank T. Ishii for pointing 
out this error. 
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from which we see that the agreement with our result (13-431) is not very good, particularly 
for ji A o. This is because the contributions of \ii=\ ~ 0(N C ) and Hi=a ~ 0{1/N C ) roughly 
cancel each other out for A and hence the 1/N C corrections cannot be neglected. To make a 
more reasonable comparison, we should compare /ij=o and \ij=\ with experimental or lattice 
results as explained above. Then, the average of the magnetic moments in fl3-45[) . 

-QiA++ + HA+ + VA° + /^-)|lattice ^ 1-27 , (346) 

should be compared with our result for the isoscalar component |/i/ =0 (Z\) = |/i/=o(p) 
~ 1.26. The isovector components are extracted by considering the differences: 

^ {HA++ - HA- ) | lattice - 248 , i (fJ, A + - fl A - ) (lattice - 2.47 , 

(HA° - Ha-) I lattice ^ 2.51 . (347) 

These values are compared with our result ^Hi=i(p) ~ 2.82. 
3.5. Axial coupling 

As explained in Ref. [14~|) . the axial coupling g A of the baryon states of/ = J = l/2is 
given by 

J d 3 x (B\s'\J^\B,s) x 2 = lg A (o- t ) s , s (T a ) I , h . (348) 
From (12-991) . the integral of the axial- vector current becomes 

d 3 x J%< = -7r 2 /tp 2 tr(ar i a- 1 r a ) I d 3 xd 2 H A . (3-49) 
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Although this is an integral of a total derivative, it does not vanish because of the term 
proportional to 1/r in (12-951) . In fact, the integral can be performed using the Gauss' 
divergence theorem as 

/ d 3 x d 2 Y 2n = [ dS- VY 2n = 4tt lim r 2 d r Y 2n = 5 0n , (3-50) 
J Js r ^°° 

where Y n (r) is defined in (I2-73P with Ao = for n = 0. This implies that only the n = 
component of the mesons, that is the pion, contributes to the integral. Again using f)3-29p . 
we obtain 



/ 



d*x (B',s>\Jr\B,s) = -^(-£-\ 



° l )s>s{T a )lL h , (3-51) 
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where (p 2 /k(Z)) is the expectation va^ 
paring this with (!3-48j) . we obtain^ 



ue with respect to the spin 1/2 baryon states. Com- 



IQttk I p 2 \ , 

°* = — {W))- (3 ' 52) 

If we approximate p and Z by their classical values, we obtain 

2iV c [6 . 
A ~ _£«/-- 0.697. 3-53 

07T V 

If we use the wavefunction f!2-22p and (12-31) to numerically evaluate the expectation value for 
nucleons, we obtain 

(m) S!lM A- (3 ' 54) 

and 

g A ~ 0.734 . (3-55) 

The experimental value 23 ' and the prediction obtained from the Skyrme model in Ref. [T4"|) 

are 

SuUp ^ 1-27 , g A \ANw - 0.61 . (3-56) 

For the excited baryons iV(1440) (Roper) with (n p ,n z ) = (1,0) and iV(1535) with 
(n p ,n z ) = (0,1), by evaluating (p 2 /k(Z)) using their wavef unctions, we obtain (see Ap- 
pendix |A~2] for details) 

^(1440)) „ ^ ^(1535)) „ Q (3 . 57) 

We can see that the axial coupling for iV(1440) is large while that for the negative parity 
baryon iV(1535) is small compared with the proton/neutron. 

There is another more direct way to perform the integral in (!3-49j) . which will be useful 
in §£T\ Using (12^951) and fl2~74j) . the integrand of (ERgjl is 

t oo \ oo 

-^j-Ef^ 2n(Z) U 3 (x-X)-J29a^ 2n (Z)Y 2n (r) . (3-58) 
^ ' n=l 2n J n=l 



*> The sign of this equation is taken to be positive such that the axial coupling is defined to be positive. 
This sign can be flipped if one exchanges the definitions of "left" and "right" chiral sectors in the two 
asymptotes z — > ±oo, i.e. the positive sign is a convention. 

**> This expression, once the classical value (k(Z)) = 1 is imposed, is equal to that obtained in Ref. . 
Furthermore, if we use (p 2 ) = p^, it agrees with gA,mag in Ref. fT7|) . 
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Following the same logic as in f l3-3p . we can show that 

oo 

J2^ 2n (Z) = MZ) , (3-59) 



where ^o(z) = ^arctanz, using f l2-40p and (!2-76p . This relation implies 

ff a « n , 2 1 



E ^. (Z ) = i (3-60) 



n=l 

and hence 



d)H A = -Y J 9a-d z ^2n{Z)Y 2n {r) . (3-61) 

n=l 

Then, we can perform the integral as 

d 3 xd]H A = / rfrrVvW)^ r = -T7^, (3-62) 
Jo ^ Tik(Z) 

using f l3-60p . From this, it is easy to rederive (13-51 p . 

3.6. Goldberg er-Treiman relation 

The axial coupling obtained in the previous subsection is related to the tcNN coupling 
QttNN by the Goldberger-Treiman relation 

9A = ■ (3 - 63) 

We can derive this relation in our context following the argument given in Ref. [Tip for the 
Skyrme model. 

It is argued in Ref. IT4"P that the pion field behaves asymptotically as 

™-lS^ v °> < 3 ' 64 > 

in the presence of a nucleon. Here, g n NN is the ttNN coupling and the expectation value 
is taken for the nucleon. In the present case, (I2-85P provides us with information on the 
asymptotic pion field. The pion field can be read from the n = component of (12-71 p 
substituted in (F2J35]). We find 

A G Z ~ n(x)Mz) + ■■■ (3-65) 

with 
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which implies 

V\XT7)/ 



mx) )^-^(j^) X -(a-r-) . (3-67) 



Comparing ()3-64p and (l3-67p . we obtain 

g-KNN _ I p 2 \ _ 1 [% _ g A 



M N " 3 \k{Z)/ 2 V k 9A f-, 



(3-68) 



where we have used (I2-38P and (13 -52 p . This is nothing but the Goldberger-Treiman relation 
(l3~63j) . 

3.7. Axial radius 

We define p A (r) as a function proportional to the expectation value of the integrand of 
(I3-62P with the normalization dr p A {r) = 1, 

Pa{t ) = (r^g a nd z i; 2n (Z)e-^\ / (^^)) • (3'69) 
The axial radius is obtained as 



(r 2 ) A ee / drr 2 p A (r) = (F A (Z)) / {k{Z)~ l ) , (3-70) 
Jo 

where 

oo 

F A {z) = 37rJ2^9 z i/; 2n (z) . (3-71) 

n=l A2 « 

From (I2-36P and (I3 59j) . this function satisfies the differential equation 

d z (k(z)F A (z)) = -37rh(z)Mz) • (3-72) 
Integrating this equation, we obtain 

k(z)F A (z) = f Q - 3tt f dz'h(z')Mz') , (3-73) 
Jo 

where fo is a constant. In order to fix f , we use the identity 

JdzF A (z) = 0, (3-74) 
which follows from (13 -71 p . Substituting (I3-73P into (l3-74p . we obtain 



/o = 3 f dz^—l dz'h(z')^ (z') - 7.82 . (3-75) 
J-oo k(z) J 



26 



Therefore, if we approximate the expectation value by its classical value, we obtain 

(r 2 ) A = (F A (Z)) I (k(Z)- 1 ) ~ F A (0) = / ~ 7.82/M| K ~ (0.582 fm) 2 . (3-76) 

If we numerically evaluate the expectation value of F A (Z) using the wavefunction (12-221) . we 
obtain 

<r 2 )^ /2 ~ 0.537 fm . (3-77) 

The experimental value is (r 2 ) 1 / 2 | cxp — 0.674 fmP^l*) 

We can compute the axial radius for excited baryons in the same manner as for the other 
quantities evaluated before. The computation of the axial radius is independent of the profile 
of R(p) in the wavefunction; thus, for the 7V(1440) (Roper) with (n p , n z ) = (1, 0), it gives the 
same result as that for the proton/neutron (13- 77p . For the iV(1535) with (n p ,n z ) = (0,1), 
we obtain (see Appendix IA.2I for details) 

(r 2 )^ 2 ~ 0.435 fm . (3-78) 

This is smaller than the axial radius of the proton/neutron. 

§4. Form factors 

In this section, we compute the form factors of spin 1/2 baryons associated with the 
currents Jy A obtained in the previous section. To this end, we first give a brief review of 
how to compute the form factors from the matrix elements of the currents. The present 
model enables us to compute the matrix elements easily with the tools formulated in §2J 

4.1. Formalism 

Let us first consider the matrix elements of a vector current for a baryon of spin 1/2, 

(p',B',s'\J^(0)\p,B,s). (4-1) 

Here \p, B, s) and \p , B', s') denote the initial and final states of the baryon under consid- 
eration. In the present section, we focus on the case where B and B' have the same n p and 
n z while the isospins /3 and F 3 may be different. The states are normalized as in f)2-2ip . The 
most general form of the matrix elements consistent with the symmetries and conservation 
of the current is 

(p',B',s'\J^(0)\p,B,s) = i(27rr 3( ^^u(p' : s f )r^(p',p)u(p,s) (4-2) 

*' This value is obtained by applying the formula (r 2 )^ = — 6d/(dk 2 ) log gA(k 2 )\k2 =0 to the axial form 
factor g A {k 2 ) in Ref. [Ml): which is fitted by a dipolc. 
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with 



^F^k 2 



2m 



u{p',s')r^ 2Z) (p\p)u(p,s) = u(p',s') 



l^F^k 2 ) 



-a^k v F 2 {k 2 ) 
a^k v F 2 (k 2 



2rriB 



u(p,s), (4-3) 
u(p,s) , (4-4) 



where k = p — p' and mg is the baryon mass. On the right-hand side of (14-21) . no summation 
is taken for the index C. Fi(k 2 ) and F 2 (k 2 ) are the scalar functions of k 2 called the Dirac 
and Pauli form factors, respectively, whose dependence on n p and n z is indicated implicitly. 
They can be computed by evaluating the matrix elements of the current. u(p , s) and u(p' , s') 
are the Dirac spinors associated with the initial and final states of the baryon of mass mg, 
respectively. The normalization conditior0 is given by 



u(p , s')u(p, s) = 5, 



m B 

-.' s n 



(4-5) 



It is useful to write the matrix elements in the Breit frame with p = —p = k/2 and 



E = E' = \lm\ + k 2 /4: 



(--,B',s'\? v (0) \^,B,s) = (2nr 3 ^5 s , s ^G E (k 2 ) , 

~ B\ s'\Jlm ~ B, s) = (2vr)- 3 ^ ±e Jla h (a a ) s , s G M (k 2 ) 



k 
2 
k 
~2 
k 
2 
k 
2 



B',s'\J${0) \'-,B,s) = (2tt 



k 
2' 
k 
2' 

k 



h r rriB 



2 
k 
2' 



2 E y 



[--, B>, s'\J c v \0) | ■-, B, s) = (27T)- 3 ^e jla h (a% s G M (k 2 ) 



(4-6) 



Here GE,M{k 2 ) are the Sachs form factors, related to the Dirac and Pauli form factors by 



G E (k 2 ) = F^k 2 ) - F 2 (k 2 ) , G M (k 2 ) = F 1 (k 2 ) + F 2 (k 2 ) , 
4m B 



G E {k 2 ) = F,{k 2 



k 2 



F 2 {k 2 ) , G M {k 2 ) = F l {k 2 ) + F 2 {k 2 



(4-7) 



The formulae of the Dirac spinor needed for this manipulation are summarized in Appendix 
IB. 21 The Sachs form factors can be obtained by evaluating the left-hand side of (14- 6p using 
the baryon wavefunctions given in §2.11 and Appendix IA.1I 



*) With this normalization, we assign \/2p° u(p, s) and \/2p° u(p, s) to an incoming and outgoing 
external line, respectively, in the computation of the Lorentz invariant matrix elements. 
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The electromagnetic form factors are defined by considering the matrix elements of the 
electromagnetic current, 



(4-8) 



Then, for the states with J 3 = +1/2 (p) and 1% = —1/2 (n), the Sachs form factors associated 
with the electromagnetic current are given by 



G P E , M (k 2 ) = 2 ( +G E , M (k 2 ) + — G EM (k 2 ' 



-G EM (k 2 ) + —G EM (k 2 )) , (for/; 



(for h = +1/2) 
-1/2) 



(4-9) 



respectively. and G v ^ are called the electric and magnetic Sachs form factors, and 

their Fourier transformation provides the distribution of the electric charge density and the 
magnetic current density, respectively. (See Refs. 1371) - f39|) for reviews.) 

Next we study the axial form factor associated with the axial current J A [x). As in the 
vector current case, we consider the matrix elements of the axial current for a spin 1/2 
baryon, 

{p',B',s'\J^(0)\p,B,s) . (4-10) 

The matrix elements consistent with the symmetries can be written in terms of the axial 
form factor #a(^ 2 ) and the induced pseudoscalar form factor gp{k 2 ) as 

-C\ 



(p',B',s'\J^(0)\p,B,s) = (2nY 



u{p',s')r^ c) (p',p)u(p, s) 



with 



u(p', s')r% (0) (p',p)u(p, s) = u(p', s') \i^^g A (k 2 ) + ^— F75 gp(k 2 )] u{p 



2m 



B 



u(p',s')r^ {123) (p,p)u(p,s) = u(p',s') hsl* 1 9A(k 2 ) + ^-k^^gpik 2 ) u(p, 



(4-11) 

(4-12) 
(4-13) 



On the right-hand side of (14 • 1 1|) . no summation is taken for the index C. The current 
conservation law yields 



gp(k 2 



g A {k 2 ), g P (k< 



Am 2 B 

k 2 aAV " 1 ' ' ~W 

In the nonrelativistic limit, the spatial component becomes 

%3 



g A {k 2 ) . 



(4-14) 



(p',B',s'\J> A (0)\P,B,s)~-(27i) 
(p',B',s'\J c A \0)\p,B,s)^-(2n) 



-3 1 _a 



kjk a 



kjk a 



Mk 2 ) , 
9A{k 2 ) 



(4-15) 
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Again, the axial form factor 5u(^ 2 ) can be computed by evaluating the matrix elements in 
(14- 151) from the baryon wavefunctions given in §2.11 and Appendix IA.1I 

In order to derive the form factors from the above formalism, it is useful to perform the 
Fourier transformation of the currents defined by 



d 3 xe- ik - 3 J"(x) . 



(4-16) 



Using the explicit form of the current flZjgg) , flMTjl . (frMty . and flZ31D , it is not difficult to 
show that 



J v \k)=e-^ Nc ^ 9 ^ 2n - l{Z) 



tli v 

2 ^ 2 + A 2 „-i 
Jv {k) ~ 6 — 1 Mo + 16*** 



2,, ^ 



i k 2 + X 2n -i 



+ 



(4-17) 
(4-18) 



Jy\k) = e 



-ik-X 



I c - i 2tx 2 k P 2 tr (T c aT a a- 1 ) e jU k x -f- 

Mq 



E 



n>l fc2 + ^2n-l 
9v^2n-l(Z) 



+ 



n>l ^ + ^2n-l 



(4-19) 
(4-20) 



= e-« (-1) J^- J a ( S n - ^i) £ g^jMg) 4 



327r 2 fi; 



A' 2 



n> 



2/i 



J^(Af) = e -*^2vr 2 /€p 2 tr (r^o" 1 ) ( <5 aj 



k a kj 
P 



E 



ga"d Z ll>2n(Z) 



+ 



(4-21) 
(4-22) 



n >1 & 2 + A 2 « 

where we have used (I2-I8p . ()2-47p . and ( l3-2ip . Here '■ • • ' denotes the terms that are odd 
with respect to they do not contribute to the result for the matrix elements. Useful 
formulas here are 

1 



d 6 xe- ik - s Y n {\x- X\) = -e 



-ik-X 



k 2 + X r 



d 3 xe- Ss H A (Z, \x-X\) = -e-** 1 ^d z ^ 2n {Z) 



(4-23) 
(4-24) 



k 2 n=l k 2 + X 2n 

The latter can be shown by using (13-60 p . In addition, it is important to note that the operator 
ordering in the Fourier transformation of the currents is fixed uniquely by the requirement 

(4-25) 



> ~ J . . • . . 

*' More precisely, Jv includes a term that is even in Z and proportional to Z. We discard this term, 

since it is negligible for large iV c and A. 
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which is equivalent to the Hermiticity of the current in the ^-representation. 

Using these expressions for the Fourier transformed currents together with the relation 
(p'\e~ lk ' x \p) = 5 3 (k — p + p'), the matrix elements can be calculated as 

/d 3 k ~ - 

— (p',B',s'\J^k)\p,B,s) . (4-26) 

4.2. Form factors 

Comparing fliTTjl - ffl^Uj) with flM}, we find 

G E (P) = N c y j 9v : { ' p2n - l{Z)) , G M (P) = N c 9 -^T 9 ^ 2n - l{Z)) , (4-27) 
„>i k 2 + X 2n .i 2 n>i A; 2 + A 2n _ 1 



G ^ = r g^Mg) Gm(^) = ^V ^ 2 "- i(Z)) , (4-28) 



with 



m s 327r^m B 2 



being the isoscalar and isovector g factors of the baryon B, as derived in § §3.21 and I3.4[ 
respectively. Here we have used the formula ()3-29p . We note also 

G E (0) = N c , G E (0) = 1 . (4-30) 

Using (14-91) . the electric and magnetic Sachs form factors for nucleons and their excited states 
with I = J = 1/2 are obtained as 

G m=E gv t 2 V (z)) » G ^ 2 )=°' 

n >l K + A 2n -1 

^) = ^E g #r M ' (4 ' 31) 

where 

» P = ^ (<7/=o + </i=0 , 9n= ^(91=0- g 1=1) (4-32) 

are the g factors of the nucleons (and their excitations). It follows that they satisfy the 
relation 

GUk 2 ) = -Glfi 2 ) = -GUP) = £ ^f 2 7 l(Z)) , GUk 2 ) = . (4-33) 
9p 9n k 2 + X 2n -i 
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Experimentally, the Sachs form factors for the proton and neutron are known to be well 
described as 

g p (?h « = V)J 1 + ^ ; Gl Ce)=0, (4-34) 

/ip fX D \ A J 

with A 2 = 0.71 GeV 2 . That is, the three form factors G P E , G P M , and G\ f are proportional 
to each other and characterized by the dipole behavior. Furthermore, the electric charge 
density of the neutron can be well approximated to be flat. It turns out that our result 
(I4-33P is in accord with these experimental results. In particular, the infinite sum in fl4-33j) 
can be approximated by a single dipole factor, showing the agreement as a function of k 2 . 
To see this, we expand our result (I4-33P as a Taylor series in k 2 , 

£ ^rf 2 T l(Z)) = (fo(Z)) - (h(Z)) k 2 + (f 2 (Z)) (k 2 ) 2 - (/ 8 (Z)> {k 2 f + ■ ■ ■ , (4-35) 

n>l k + A 2n-1 

where the coefficients are obtained as the expectation values of 

MZ) = S (A 2 „_0*« • (4 ' 36) 

n>l 

As seen in (13-31) . we have (f (Z)) = 1. Let us evaluate these coefficients using the classical 
approximation (ip2n-i(Z)) ~ V ; 2n-i(0). The first nontrivial coefficient, /i(0), is in fact equal 
to Fq/G given in f l3-10p . The higher coefficients can be obtained in the same way as F , as 
they satisfy the recursive relation 

- d z (k(z)d z f k (z)) = h{z)f k ^{z) , f k (z) = f k {-z) , / fc (±oo) = . (4-37) 

We thus obtain 

/i(0) = 2.38 , / 2 (0) = 4.02 , / 3 (0) = 6.20 , / 4 (0) = 9.35 , / 5 (0) = 14.0 , • ■ ■ . (4-38) 
On the other hand, the dipole expression ( I4-34|) can be expanded using A 2 = 0.758 GeV 2 as 

1 + ^\ = 1 - 2.38k 2 + 4.24(£ 2 ) 2 - 6.71(£ 2 ) 3 + 9.97(fc 2 ) 4 - 14.2(fc 2 ) 5 + • ■ ■ , (4-39) 

where we have chosen the parameter A 2 in such a way that the first coefficient reproduces the 
/i(0) of (I4-38P with the unit M KK = 949 MeV = 1. We find a good agreement for the latter 
coefficients, suggesting that our form factors (I4-33P exhibit the dipole behavior indicated by 
experiments. It is also useful to note that the relation (I4-37P implies that the function 

n >i k + A 2 „-i 
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satisfies 



d z (k(z)d z F(k 2 , z)) = k 2 h(z)F(k 2 , z) , F(k 2 , z) = F(k 2 , -z) , F(k 2 , ±00) = 1 . (4-41) 
The solution of this equation evaluated at z = is plotted in Fig. [TJ 




0.5 



1 . 5 



k 2 [GeV 2 ] 

Fig. 1. A plot of F(fc 2 ,0), which is equal to G^(fc), with (ip2n-i(Z)) approximated by the classical value 
V'2n-i(0). Our result (solid line) reproduces the dipole behavior (dotted line) in (|4-39[) with A 2 = 
0.758 GeV 2 . 

The electric and magnetic charge radii can be computed from the first coefficient of the 
form factor expanded in powers of k 2 . Namely, they are given by 



d 



( t2 )e.m = -6-=r-logG EM (k 



E ' M dk 2 

except for the neutron charge radius, which is defined by 



(4-42) 



fc 2 =0 



/£ ' n dk 2 M 



(4-43) 



fc 2 =0 



Since all the form factors are proportional to each other (except for G E , which vanishes) 
as in (14- 331) . we conclude that 



(r 2 ) M , P = (r 2 ) M , a = (r 2 ) E , p , <r% n = . 



(4-44) 



It is easy to check that the definition of the charge radii in (I4-42|) is consistent with our 
previous calculation in §31 and using our result for the electric charge radius of the proton 
(13-121) . we obtain 



2\V2 
■P 



„2\V2 



' M,n 



0.742 fm . 



(4-45) 



The values observed in experiments are (r 2 )^ 2 | exp — 0.875 fmpS (r 2 )]^ 2 
0.855 fm,S3 and (r 2 ) 1 ^ | cxp — 0.873 fmj^ 1 which are reasonably close to our result. 



M,p l ex P 
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The axial form factors are given by 

g A $ 2 ) = ^V M^ , ^(fc 2 ) = ^(p 2 )V ^^ 2 " (Z)) . (446) 
Note that the values at k 2 = are obtained by using (I3-60P as 

~ , s N c / 1 \ , . 167TK / p 2 \ , 

^ (0) =^(i(I))' ^ (0)= — (ife)- (447) 



16vr 

The latter reproduces (13- 52() . 

It is empirically known that the axial form factor gA{k 2 ) can also be well fitted by a 
dipole profile. Using the same technique as above, the Taylor expanded axial form factor 
with the classical approximation (dzip2n(Z)) ~ d z ip2n(0) is found to take the form 

~ 1 - 1.30fc 2 + 1.09(£ 2 ) 2 - 0.770(£ 2 ) 3 + 0.511(£ 2 ) 4 - 0.331(P) 5 + ■ • • . (4-48) 

9a{0) 

These coefficients are close to those obtained from the dipole profile 
1 + P 



M\ 



= 1 - 1.13A; 2 + 0.958O 2 ) 2 - 0.721(A; 2 ) 3 + 0.510(£; 2 ) 4 - 0.345(fc 2 ) 5 + • • ■ , (4-49) 

with Ma 1.26 GeV. 
4.3. Cubic coupling 

The form factors computed above are composed of an infinite tower of poles that cor- 
respond to the vector and axial-vector meson exchange. From the residues, we can extract 
information on cubic couplings among baryons and (axial-)vector mesons. 

We assume that there exist cubic couplings of the form 

n>l ^ ' 

+ 4^~ J2 (J^bb (dffi - d v %) Ba^ ^-B + h v n BB (8X a ~ duv; a ) Bo*" y£) (4-50) 

n>l ^ ' 

in the four- dimensional baryon-meson effective action. Here v^(x) and v™ a (x) are the U(l) 
and 57/(2) parts of the nth vector meson v™(x), respectively, and B(x) is the baryon field. 
Note that we do not include the direct interaction among baryons and the background gauge 
potential Vj^~\ since the present model exhibits the complete vector meson dominance as 
argued in Ref. [6]) and emphasized in Ref. [17]) for the cases including baryons. 
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We recall that the vector and axial- vector mesons couple with the background gauge po- 
tentials as in f )2-34p . Using these interactions, the Dirac and Pauli form factors are computed 
to be 

£ ( £ )= £*^™ £ ( £ ) = £*^» , (4 .51) 

n > x K + A 2 r»-1 n >x « + A 2n -1 

and a similar relation holds for the SU(2) sector. On the other hand, in the present model, 
(I4-27P and ( I4-28P together with (14- 7p yield the Dirac and Pauli form factors of the form 

F 1 Ck 2 ) = N c y2 9 ^ 2n - l{z)) , fS 2 ) = nJ 9 -^-i)J2 sA^{z)) ? 

n >i A; 2 + A 2n _i ^2 J n >! k 2 + A 2 „_i 

F 1 (fc 2 ) = E ^ 2 "' l(Z)) . F 2 (k*) = 9 ^T 9 ^ 2n - l{Z)) . (4-52) 

n >! fc 2 + A 2rt _! 2 ^ fc 2 + A2n-l 

Here we have kept only the leading terms in the large N c and large A limits 
By comparing these results, we obtain 

dv^BB = N c (lp 2n -l{Z)) , h v n BB = N c ~ lj (fan-l {Z)) , 

9v n BB = {lp2n-l(Z)} , h v n BB = ^yi (ll) 2 n-l(Z)} . (4-53) 

It is interesting to note the relation 

g V "BB = N c g v n BB . (4-54) 

In particular, for the case with n = 1 and B being a nucleon, this means the g^NN = N c g P NN, 
which coincides with the constituent quark model prediction. This relation is found also in 
Ref.dZ]). 

The axial form factors ^4 and g& can be obtained from the effective cubic couplings with 
pion and axial-vector mesons of the form 

Ct = E (g^BB^Bi lbl ^B + g a n BB a™ Bi 7 57"yfl) 

n>l ^ ' 
( — T° — T a \ 

+ 2i ( g^BB 7T Bj 5 —B + g nBB 7r a Bj 5 —B j , (4-55) 

where a^(x) and a™ a (x) are the U(l) and SU(2) parts of the nth axial vector meson a^(x), 
respectively, and tt(x) and ir a (x) are the U(l) and SU(2) parts of the pion field II(x), 



*' Since the leading contribution to the baryon mass tub is Mq = 8it 2 k, we consider tub to be of order 
AiV c . However, we will not use the relation mj ~ Mo, since we know that the subleading contributions in 
m b are not small as discussed in Ref . [TB| . 
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respectively. As in (14-501) . no direct coupling between the baryons and is assumed. It 
follows from this and (!2-34[) that the invariant amplitude with the incoming and outgoing 
baryons in the presence of the external gauge field V/i ^ is given by*3 



^757/, J2 



9a n 9a n BB 

k 2 + X 2n 



n>l 
2f- K g 1 , BB 



k 2 



9a"9a n BB 



n>l 



A 



2// 



k 2 + X 



2/i 



u(p,s) 



+ (SU(2) part) . 



(4-56) 



By comparing this with the matrix elements (14 -111) (14 -131) . we obtain g^, h^, 9a-, an d as 
functions of the cubic coupling constants in (I4-55|) : 



9A(k 2 ) = y: 



n>l 



9a n 9a n BB 

k 2 + A 2n 



9 a 



(* 2 ) = £ 



n>l 



9a n 9a n BB 

k 2 + X 2n 



(4-57) 



g P {k 2 ) = 2m B 



g P {k 2 ) = 2m B 



2/^ 



BB 



k 2 

2 fn9iTBB 



4m 2 B J2 



9a n 9a n BB 1 



n>l 



A 



2ii 



k 2 + X 



2/i 



k 2 



4m 2 B J2 



9a n 9a n BB 



n>l 



A 



2n 



k 2 + X 



2// 



(4-58) 



Equating these results with ( I4-46P then using ( I4-I4p and (I3-60I) leads to the relations 



9a n BB 



32lT 2 K 



(dzip2n(Z)) , g a * 



^ 2 K 



BB 



(p')(d z ip 2n (Z)} 



(4-59) 



^ m B N c I 1 \ m B 167TK / p 2 \ 

9nBB = J^YO^ \k{Z)/ ' 9nBB = 

Comparing (I4-6U|) with (I4-47|) . we note that the following relations hold: 



9a(0) 



fir9TT 



BB 



m B 



9a(0) 



fir9TT 



BB 



m B 



i.e., the Goldberger-Treiman relation. 
4.4. Numerical estimate 



(4-60) 



(4-61) 



By solving ( I2-36P numerically using the shooting method, the expectation values of 
ip2n-i(Z) and dzip 2n {.Z) with respect to the wavefunction ipz(Z) (see Appendix I A. 1 1 for 



The propagator of a™ is given by that of a Proca field: k -i+ x (Vnv + fyu&i>/A2n)- 
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details) can be estimated as*3 



n 


An 


(ip n {Z)) nz=0 


{lp n {Z)) nz=l 


(dz4>n{Z))n z =Q 


(d Z ^n(Z)) nz = 1 


1 


0.669 


5.80 


4.51 








2 


1.57 








3.46 


0.618 


3 


2.87 


-2.70 


0.766 








4 


4.54 








-3.08 


2.22 



(4-62) 



As an example, consider the nucleon with n. 



n 7 



0. Using this table, we obtain 











n 


9v n NN 


Qv n NN 


9a n NN ga n NN 


1 


17.4 


5.80 


4.42 6.14 


2 


-8.10 


-2.70 


-3.84 -5.46 



(4-63) 



In particular, for the p meson we thus obtain 

QpNN = 9v 1 NN 



5.80 



(4-64) 



This is consistent with the experimental data g P NN\cx P = 4.2 — 6.5® (see also Ref. ITTj) ). 
Note that the functional form of (I4-53P shows the universality of the p meson couplings 
among the spin 1/2 baryons with any n p and n z = 0. Moreover, for large N c and large A, 
9v 1 bb tends to take a common value for any baryon state B with / = J = 1/2, because 



the corresponding wavefunction has a narrow support with a width of 0(X~ 1 ^ 2 N c ). For 
reference, 5^(1535)^(1535), the cubic coupling of p with iV(1535), can be numerically computed 
using the quantum wavefunction for n z = 1 as 



67^(1535)^(1535) — 



n, = l 



4.51 



(4-65) 



To see if the p meson universality holds extensively in the meson sector, we list the meson 
cubic couplings g p7r7r , g P v n v n , and g pa n a n computed in Ref. [6j): 



gpTTTT 



4.81 



n 


g P v n v n 


g P a n a n 


1 


5.19 


3.32 


2 


3.12 


2.98 


3 


2.93 


2.89 


4 


2.87 


2.85 



(4-66) 



It seems that the relation g p7T7T ~ g ppp ~ g P NN roughly holds within 20% error, although 
g pv n v n with n > 1 and g pa ™ a n are not close enough to ensure the universality. 

*' Here we show values only up to n = 4, as the quantum treatment of the values seems to break down 
for higher n because its deviation from the classical value becomes significant. 
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If we use tub = 940 MeV as an input, together with f l2-3|) . the Yukawa couplings of the 
pion and the nucleons given in (l4-60p are evaluated as^ 



9nNN ^ 5.37 , Q-kNN ^ 7.46 



(4-67) 



while the experimental value is gvAnvlexp — 13. 2.^ The smallness of the predicted value 
is related to the similar observation for the axial coupling analyzed in §3.51 through the 
Goldberger-Treiman relation (14-611) . In the same manner, the Yukawa couplings involving 
iV(1440) and iV(1535) are easy to evaluate: 



5 l 7rAf(1440)A r (1440) — 8.23 , 5^(1440)^(1440) — 16-7 , 
SViV(1535)jV(1535) — 4.55 , S , 7 rJV(1535)Ar(1535) — 6.32 . 

Here we have used fTAd3l) . flA44j) . flA46]) . and (lAd7j) . 



(4-68) 



**) 



§5. Summary and discussion 

In this paper, using the model of holographic QCD proposed in Refs. E]) and E]), we 
calculated the static quantities of the nucleon and excited baryons such as iV(1440), iV(1535), 
and A, and the form factors of the spin 1/2 baryons. The baryons are described as quantized 
instantons in five-dimensional YM-CS theory.^ By defining the chiral currents properly at 
the spatial infinity in the fifth dimension and by solving the YM-CS equations of motion 
using the instanton profile given in Ref. [IB"]) and Green's functions, we obtain an explicit 
expression for the chiral currents depending on the baryon state. From the currents we 
computed various static quantities of the proton/neutron, such as the charge radii, magnetic 
moments, axial coupling and axial radius. See the summary table below. The Goldberger- 
Treiman relation is naturally derived. These quantities can be computed for excited baryons 
in the same manner, which are our theoretical prediction for the excited baryons (see the 
second table below). We also calculated the nucleon form factors (and also those for excited 
baryons). It was shown that the electric and magnetic form factors of the nucleon are 
roughly consistent with the dipole behavior observed in experiments. The electric as well as 
the magnetic charge radii of the baryons and their couplings to mesons are calculated from 
the form factors. 

*' drrNN is related to the Yukawa coupling of the isoscalar mesons g^NN and g^'NN ■ It is, however, 
difficult to tell which component it corresponds to, since we are analyzing massless QCD with Nf = 2. 

**> A computation of the Yukawa couplings involving the excited baryons is also performed in Ref. 
where a five-dimensional spinor field is incorporated in the gravity side in the bottom-up approach in holo- 
graphic QCD. 
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A table summarizing our results for the static properties for the proton and neutron is 
given below. For a comparison, the table includes the values obtained in experiments, and 
also the results obtained from the Skyrmion.^ 





our model 


Skyrmiorf® 


experiment 


V /M,I=0 

( r2 )M,p 
( r2 )M,n 

(r 2 )f 


0.742 fm 
0.742 fm 
(0.742 fm) 2 


(0.742 fm) 2 
(0.742 fm) 2 
0.537 fm 


0.59 fm 
0.92 fm 

OG 

— oo 

OG 
OO 


0.806 fm 
0.814 fm 
(0.875 fm) 2 
-0.116 fm 2 
(0.855 fm) 2 
(0.873 fm) 2 
0.674 fm 








2.18 


1.87 


2.79 




fin 




-1.34 


-1.31 


-1.91 








1.63 


1.43 


1.46 




9 A 




0.734 


0.61 


1.27 


Q-kNN 


7.46 


8.9 


13.2 


QpNN 


5.80 




4.2 ~ 6.5 



For excited baryons, iV(1440) (Roper) and iV(1535), we provide our theoretical predic- 
tions in the following table. 
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n,p 


iV(1440) 


iV(1535) 


(r 2 ) E 


>p 


(0.742 fm) 2 


(0.742 fm) 2 


(0.699 fm) 2 


(r 2 ) E ,n 











(r 2 ) M , P 


(0.742 im) 


/n ^7 jo r \2 

(0.742 im) 


(0.699 tm) z 


( r /M,n 


[U. 1 4/ M ) 


('PI 7/10 fm"\ 2 

l^u. < 4/ im ) 


^u.oyy im j 




U.OOI III1 


n ^7 fm 

u.ooi im 


n /I ^ fm 

u.4oo im 




flp 




2.18 


2 99 


2.18 








-1.34 


-2.15 


-1.34 




Pn 




1.63 


1.39 


1.63 




9a 




0.734 


1.07 


0.380 


9-kBB 


7.46 


16.7 


6.32 


9pBB 


5.80 


5.80 


4.51 



As shown in the first table, we found good agreement with experiments for various 
quantities of baryons. Our numerical results presented here should be treated with caution 
since they were obtained for large values of the 't Hooft coupling and N c . In order to 
incorporate the difference among excited baryon states, we included subleading corrections 
only at the last stage of the calculations. This procedure is difficult to justify since we 
considered the leading order action ( 12- ip as our starting point, and there should be more 
corrections. Further investigation to improve the accuracy of the results would be interesting. 

In this paper, baryons have been described as a soliton, and the baryon physics has been 
analyzed by a standard semiclassical quantization of the soliton. This approach appears to 
be completely different from the treatment of Ref. ITTp . in which a key step is to introduce 
a five- dimensional spinor field into the five-dimensional YM-CS system to represent the 
baryons. In order to relate the two approaches, consider the vBB coupling, for simplicity, 
in the four-dimensional meson-baryon effective Lagrangian (I4-50|) : 



hit 



£ 

n>l 



g v n BB v^{x) B(x)i^ 1 —B{x) + 



(5-1) 



If we substitute g v ^BB from (14-531) . this can be written as 



int 2 



dz [A IJ ,(x,z)B(x,z)ij tJ '—B(x,z) + 



(5-2) 



Here is the U(l) part of the five-dimensional gauge potential in f )2-35p . and B(x, z) = 
B(x)ipz(z) is a five-dimensional spinor field constructed by the four-dimensional baryon field 
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B(x) and the suitably normalized wavefunction ipz{Z) corresponding to the baryon B. This 
term can be summarized as a five-dimensional gauge interaction by regarding ipz{Z) as 
an eigenmode of a wave equation that reproduces the baryon spectrum found in Ref. [I~6j) . 
Likewise, the interaction terms given in (14-501) and (I4-55P are sufficient to reconstruct the 
five-dimensional spinor field action with the Pauli-type interaction in Ref. IPT|) . 

We have concentrated on a one-point function of mesons and the electromagnetic field in 
the presence of a single quantized soliton, mainly to extract the static properties of baryons. 
There are other interesting aspects of the force associated with the baryons, for example, 
interactions between baryons, in particular. When baryons are far from each other, we can 
use our results and compute the nuclear force based on a one-meson-exchange picture of the 
nuclear force. Resolving the issue of the nuclear force at short distance is interesting and 
will be reported in our forthcoming paperp3> 
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Appendix A 

Excited Baryons 

In this appendix, first we summarize quantum wavefunctions for the lowest /excited states 
of baryons obtained in Ref. [T6l) . and present detailed calculations relevant to the static 
quantities of the baryons presented in 

A.l. Wavefunctions for excited baryons 

In Ref. [T6il . fluctuations around the classical solution (I2-I2p of the five- dimensional YM- 
CS theory in the curved background (12- ip are quantized. The quantization of the soliton 
was done in the approximation of slowly moving (pseudo)moduli of the soliton solution. The 
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moduli degrees of freedom of the soliton are 



X\t) , Z(t) , p(t) , a\t) . (A-l) 

X 1 and Z describe the center-of-mass motion of the soliton, while p and a 1 (I = 1,2,3,4) 
describe the size of the soliton, which is instanton-like in the SU(2) sector, and the orientation 
of the instanton in the SU(2) group space, respectively, with (a 1 ) 2 = 1. The Hamiltonian for 
the moduli (I2-17P provides the quantized energy eigenstates of the baryon, specified by the 
quantum number B = (/, I 3 , n p , n z ) and its spin s. The wavefunctions for the quantized states 
can be written explicitly for low-lying states. The nucleon wavefunctions, B = (1, ±1/2, 0, 0) 
and s = 1/2, are written as ( 12 -22 p . 



pt) ^ R(p)ipz(Z)(a 1 +ia 2 ) , \nt)^R(p)^z(Z)(a 4 + ia 3 ), (A-2) 



with 



R(p) = p-^V^^ e -^ p2 , MZ) = e~^ z2 . (A-3) 

Here we present the first excited states for which static quantities are computed in £j3j 
The excited state with B = (1, ±1/2, 1, 0) corresponds to iV(1440), which is called the Roper 
excitation. The wavefunction is 



R{p) = (^p 2 - 1 - 2Vl + JV?/5) p-^^^e-^ 2 . 



(A-4) 



The excited state with B = (1, ±1/2, 0,1) is iV(1535). The wavefunction for the Z part is 
now given as 

MZ) = Ze-^ z2 , (A-5) 

where, again, its normalization constant is not fixed yet. The other part of the wavefunction 
is the same as that of the proton/neutron. This excitation has a negative parity, which is 
reflected in the oddness of the wavefunction ip z (Z). 

The excited states in the SU{2) group space have already been studied in the context 
of the Skyrmion. 14 ' The lightest among these excited states is A with I = J = 3/2. The 
isoquartet is composed of the four baryons A + , A , and A~ with J 3 = 3/2, 1/2, —1/2, 
and —3/2, respectively. As shown in Ref. [T4"]) . A ++ with s = 3/2 is described by the 
wavefunction 

(ai + ia 2 f . (A-6) 
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The rest of the A wavef unctions with s = 3/2 are obtained by letting the isospin lowering 
operator act on ( 1A-6|) : 

(ai + za 2 ) 2 (a 4 + ia 3 ) , (for A + ) 
(ai + ia 2 )(a 4 + ia 3 ) 2 , (for A ) 

(a 4 + za 3 ) 3 . (for ZT) (A-7) 
The wavefunction R(p) for A (states with I = 3 and n p = 0) is given as 

R(p) = p -l+2v/4+^/5 e -^P 2 (A . g) 

We use these wavefunctions to compute the static quantities of excited baryons. 

A. 2. Expectation values relevant to iV(1440) and iV(1535) 

First, we summarize what quantities are necessary for computing various static quantities 
presented in $3) 

(F(Z)) , (p 2 ) , (j^) , (F A (Z)) . (A-9) 

The first quantity F(Z) is necessary for computing (r 2 ) I=0 , and thus (t 2 )e,m- The second, 
third, and fourth quantities, are for the magnetic moment /x, the axial coupling g A , and the 
axial radius (r 2 ) 1 / 2 , respectively. 

First, we compute (p 2 ). The integral to be evaluated is 

{p ^ = JpMpNp ' ( } 

and we substitute (1A-4P into this to allow for the Roper excitation. By partial integration, 
we obtain 



2 y/5 + 2 v /5TiVf 

(p )n p =0 = Pel ' ^ ' 

2 3y^+2 V / 5TiVf 

(p )n p =l = Pel • ( A ' 12 ) 

The former equation is for the proton/neutron, (13 -33 p . For N c = 3, these are numerically 
given by 

(p 2 ) np=0 ~ 1.62 x p 2 cl , (A-13) 

(p 2 ) np=1 ~ 2.37 x p 2 cl . (A-14) 
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Next, we compute \^z)J^ (F(Z)), an d {}a{Z)), that is, the quantities relevant to the Z 
directions. This is necessary for the iV(1535) excitation. The classical values are given with 
Z = as 

(jP))- 1 ' (F(Z)>~F(0), (F A (Z)) ~ F A (0) . (A-15) 

In the large A expansion, the wavefunction for Z is localized at the origin Z = 0. Inclusion 
of the subleading terms causes differences in baryon states. 

Numerical evaluation using the wavefunction (IA-5j) is performed as follows: 



V J ~ z " -~ 0.649, (A-16) 



/ 1 \ fdZ^e^^ 



Wz)/ nz=1 j d zz^ z2 



0.337 , (A-17) 



where we have used k = 0.00745. 

For F{Z) and F A (Z), we first solve the differential equations satisfied by them, (I3-7|) 
and 03-72 p . and then use the solutions to numerically evaluate the normalized integrals. We 
obtain 

(F(Z)) nz=0 ~ 12.7 , (F(Z)) nz=1 ~ 10.9 , (A-18) 
(F A (Z)) n ^ =0 ~ 6.67 , (F a (Z)> Bb=1 ~ 4.38 . (A- 19) 

Classical values are F(0) = 14.3 and Fa{0) = 7.82; thus, the ratios to the classical values 
are given as 

^f(o;T =0 ~ 0,892 ' ^ F ^fJo) z=1 ~ 0,762 ' ^ A ' 20 ^ 

{FA l Z h= * 0-852 , {FAi f] )n ' =1 * 0.560 . (A-21) 

Thus, smaller values of the quantities are obtained for n z — 1. 

Using these expectation values, the static properties of A^(1440) and AT(1535) can be 
computed, as presented in $3j 

A. 3. Magnetic moment of A 

Evaluation of the isoscalar and isovector magnetic moments requires the matrix elements 

of 

pY = -ip 2 tr (Va^d) , tr (a^a^V) , p 2 . (A-22) 
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Here, we compute these for the wavefunctions of A given in (1A-6|) . (1A-7I) and (1A-8j) . 

First of all, the third component of the spin of these states is chosen to be s = +3/2. 
Therefore, using (!2-47j) . we can see that the value of p 2 x* =3 is three times that of the pro- 
ton/neutron for which s = 1/2. 

Next, let us evaluate tr (ar*a _1 r a ). From the symmetric nature of this quantity, it is 
enough to consider the index i = a = 3. For this choice of the index, we have 



tr (a^a"V a ) = Aa\ + Aa\ - 2 . 



(A-23) 



In order to calculate the expectation value of this quantity with the wavefunctions flA-6j) and 
flA-7p . we use the spherical coordinates 



a 4 = cos^q , 
a 3 = sin 6*o cos 9i , 
a 2 = sin 9 sin 9i cos 9 2 
a± = sin 6q sin 9\ sin 82 



(A-24) 
(A-25) 
(A-26) 
(A-27) 



The Jacobian for this change of variables is df2 3 = sin 2 9 sin 9id9 d9id9 2 - Then we obtain, 
for instance, 

{A ++ \(Aal + Aal-2)\A ++ ) 

2tt £ £ d9 d9± sin 8 9 sin 7 9 1 (4 cos 2 9 + 4 sin 2 9 cos 2 9 l - 2) 



/ dn 3 (al + 



6 

'5 ' 



(A-28) 



{A+\{Aal + Aal-2)\A+) 

2tt £ £ d9 d9 1 sin 6 # sin 5 fli(cos 2 fl + sin 2 #o cos 2 #i)(4 cos 2 fl +4 sin 2 fl cos 2 #i -2) 



/ df2 3 (al + a\) 2 [a\ + a 



2 
5 



(A-29) 



The denominators originate from the normalization of the states (1A-6|) and ( 1A-7j) . It is also 
easy to see that 

(A°\(4al + 4a 2 - 2)\A°) = ~ , (ZT|(4a 2 + 4a 2 - 2)\A~) = jj . (A-30) 



The expectation value of p 2 with respect to the wavefunction flA-8j) is 



, 2N _ V5 + 2^WTN! 2 

x/ 5 /l=3 ~ osr /'el 



2AT, 



(A-31) 
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and the ratio to (13-331) is 



(P 2 ) 



1=3 



V5 + 2^20TN! ^ LM 



(A-32) 



for N c = 3. 

Comparing these with the proton state with up spin, we obtain ()3-4ip and ()3-42p . These 
ratios are used to compute the magnetic moments of A as 





1 . 

^ = 3 X - {fjLp + fl n 


. 9c 1 . 
) + j X 2^ 


- fl n ) ~ 5.50 , 


(A-33) 


fi A + 


1 . 

= 3x-^ + fi n ) 


3c 1 , 
+ 5 X 2^ 


- /i n ) ~ 2.67 , 


(A-34) 


fiA° 


= 3 x i (fi p + fi n ) 


3c 1 , 


fi n ) ~ -0.15 , 


(A-35) 




1 . 

= 3 x - (fip + ji n ) 


9c 1 , 


- /i n ) ~ -2.97 , 


(A-36) 



i.e., f )3-43p . Here we have substituted the magnetic moment for the proton and the neutron, 

dS3ZD. 

Appendix B 

Useful Formulae and Notation 



B.l. Useful formulae for currents 

Here we summarize useful formulae that are used in §2] for computing the currents. For 
g defined in (I2-I3p . the following equations are obtained, 

(B-l) 
(B-2) 
(B-3) 
(B-4) 

In particular, its relation to SU(2) generators is given as 

1 



gdig' 1 = 




- Z)t L - E^X 3 


- x j y) 


gdzg' 1 = 




- X b )r b , 




9~ X di9 = 




-Zy + e ijb (x j 


- X J )r b ) 


g^d z g = 




- X b )r b . 





g^g- 1 



e 



z - Z) 2 -\x- X\ 2 j r a + 
+ 2 (e bac (x b - X b ){z -Z) + (x a - X a ){x c - X c )) 



[ g- x d ig , r a ] = — ({x b - X b )5 ia - (x a - X a )c? b + (z - Z)e iab ) r b 



(B-5) 
(B-6) 
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[g-'dzg, r a ] = -h x * - X k )e kab r b , (B-7) 

gr b + r b g = ^ {(z - Z)r b - i(x b - X b )) , (B-8) 

gr b - r b g = ^{x a - X a ) e abc r c . (B-9) 
A useful identity for the epsilon tensor is 

Qjral^'l £abj% % £abi%" "I - ^ijb% % ■ (BTO) 

B.2. Useful formulae for Dirac spinors 

We summarize the formulae of the gamma matrix and the Dirac spinor, which are used 
in the previous sections. 

The Dirac representation of the gamma matrix is taken as 

7 « = -i i j °), y = ( _r * ) , T5 = iT »y 7 v =(; i 



We define 



It can be verified that 



(B-ll) 



^ = i[y*,7"] • (B-12) 



7 o r t 7 o = +7 , ; 7 o 7 t 7 o = +75; 7 o^ 7 o = _^ _ (B . 13) 
The on-shell Dirac spinor is given by 



1 VBTmix,,, , 

V2-E V yJE -m B n- ax( s ] 



with 



n = pr, (B-15) 
\P\ 

X(i/2) = J j ' X(-V2) = ] • ( B ' 16 ) 
This satisfies the Dirac equation 

m B )u(p, s) = , s)(itf+ m B ) = , (B-17) 
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with 



u = u*p. (/9 = i7°) (B-18) 

The Dirac spinor is normalized as 

u(p, s')u(p, s) = 5 ss i . (B-19) 
In the nonrelativistic limit, the Dirac spinor reduces to 

u(p,s)= ( 1 {S \ ) +0(m-/) . (B-20) 



2mg 



It is easy to verify 



u(p', s') 75 u(p, s) = -^—k a (a a ) s , s + 0{mf) , (B-21) 

ZTfl b 

HP, s') 757° u(p, s) = -?—(p + p') a (a a ) s , s + 0{mf) , (B-22) 
u(p', s') 75 y M (p , s) = i(<T j ) s , s + 0(m B 2 ) . (B-23) 

Here k = p — p' and we used 

x\ s , ) v a X(s) = (v a ) s 's. (B-24) 
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